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We observe that the fully nonlinear evolution equations of Rosenau and Hymann, often abbrevi- 
ated as K(n, m) equations, can be reduced to Hamiltonian form only on a zero-energy hypersurface 
belonging to some potential function associated with the equations. We treat the resulting Hamil- 
tonian equations by the dynamical systems theory and present a phase-space analysis of their stable 
points. The results of our study demonstrate that the equations can, in general, support both 
compacton and soliton solutions. For the K(2, 2) and K(3, 3) cases one type of solutions can be 
obtained from the other by continuously varying a parameter of the equations. This is not true for 
the K(3, 2) equation for which the parameter can take only negative values. The K{2, 3) equation 
does not have any stable point and, in the language of mechanics, represents a particle moving with 
constant acceleration. 
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1. Introduction 

The theories of dynamical systems via phase-space analy- 
sis provide a very useful tool to extract information about 
any physical system modelled by linear or nonlinear dif- 
ferential equations. Here the long term behaviour of the 
system is investigated in terms of the stability of fixed 
points of the differential equations. In particular, one 
examines how the stability or instability is affected as 
the parameter of the model or differential equation is var- 
ied. This viewpoint was followed by Caffey [1] to identify 
the closed form solution of a number of quasilinear dif- 
ferential equations with special attention to phase space 
trajectories that connect a critical point to itself or other 
critical points. 

A nonlinear evolution equation is often called quasilin- 
ear when the dispersive term of it is linear. In addition to 
quasilinear equations there are fully nonlinear evolution 
(FNE) equations characterized by nonlinear dispersive 
terms. The FNE equations play a role in the studies of 
pattern formation and wave breaking phenomena [2] . As 
a straightforward generalization of the quasilinear KdV 
equation Rosenau and Hyman [3] introduced a class of 
FNE equations given by 
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u t - A{u n ) x + {u m ) xxx = 



(1) 



Equations in (1) are often called K(n, m) equations with 
n > , 1 < to < 3. Here A is a parameter of the model. 
The FNE equations in (1) are non integrable and support 
compacton solutions. The compacton represents a trav- 
elling solitary wave solution with compact support. The 
object of the present work is to provide a phase-plane 
analysis of (1) and thus illustrate the type of insight one 
can gain from the study. 

The quasilinear equations can easily be treated by the 
dynamical systems theory. For example, the change of 
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l(x, t) = (j){x — vt) = 4>{£) 



(2) 



followed by an integration over £ converts the quasilinear 
KdV equation 
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with parameters a and 7, to an ordinary nonlinear dif- 
ferential equation 
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such that <j> plays the role of space coordinate and £, that 
of a time coordinate. Thus for stability analysis, (4) can 
be viewed as an equation of anharmonic oscillator in the 
variables <j> and £. 

2. Analogue of (4) for FNE equations 

For the FNE equations in (1), a similar change of variable 
from u to 4> and subsequent integration over £ lead to 
ordinary differential equation 
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The last term in the right side of (5) does not permit 
one to regard these equations as Hamiltonian systems 
and give rise to an awkward analytical constraint for the 
application of dynamical systems theory. This problem 
can, however, be resolved by considering evolution of <f> in 
a zero-energy hypersurface in the phase space belonging 
to some potential function. 

The potential representation for the analysis of trav- 
elling wave solutions of nonlinear dispersive evolution 
equations was introduced by Eichmann, Ludu and 
Draayer [4] . This representation is defined by 
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The left side of (6) was identified with the non-relativistic 
kinetic energy and right side with the negative value of 
a potential energy F(4>) because £ and cf) in (2) could be 
regarded as time and space coordinate respectively. Here 
^| is the velocity of a particle moving along the (/>-axis. 
Clearly, for the potential representation in (6), the evo- 
lution of 4> or u proceeds on the zero-energy hypersurface 
in the phase space belonging to The function J 7 ^) 

for the K (n, m) equation in (1) is given by [4] 
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Using (7) in (5) we get 
d 2 c 
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lead to (8). Formally, the critical points of (11) are given 

by 



where 



and 
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3. Phase-space analysis of K(n, to) equations 
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Equation (8) can be written in the newtonian form 
d 2 4> dV(<t>) 
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with the potential 
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The second-order differential equations in (8) are equiv- 
alent to two first-order equations given by 



^ = iP = V(<f>, V), (say) 
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= Q(</>, V), (say) 



(116) 

Here f/' stands for the velocity of the particle at a point 
</>(£) on the hypersurface defined by (6). It is rather 
straightforward to show that (11a) and (116) form a 
Hamiltonian system with the Hamiltonian 
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Clearly, the canonical equations 
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The travelling wave solutions of (1) are usually studied by 
restricting the values of both integers, n and to to 2 and 
3. We shall, therefore, present the phase-space analysis 
of the K (n, m) equations given by 
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with a view to study the stability of their critical points. 
The equation of motion (16d) can be viewed as that of an 
inverted harmonic oscillator while (16a) represents a sim- 
ilar oscillator with shifted origin. In the context of Bose- 
Einstein condensates the potentials in (16a) and (16d) 
are often referred to as expulsive potential [5]. Equation 
(166) gives the motion of a particle in a linear poten- 
tial. For small <f> 1 the potential in (16c) is similar to that 
of (166). For large </>, however, the potential is highly 
nonlinear. 

From (11) it is easy to see that, as with linear harmonic 
oscillator, (0, 0) is also the critical point of (16d). It is 
well known that the critical point of linear oscillator is 
a center and is neutrally stable. We shall now study 
the stability of critical point of (16d). Since all K(n, m) 
equations can be written in the form of (11) we make use 
of the eigenvalues of the Jacobian matrix 
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FIG. 1: Phase diagram for the K(3, 3) equation 



for the classification of their critical points. In (17) 
_&P u- dV \ 



and 
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where ((/>*, %/}*) stands for the critical point. For (16d) 
the eigenvalue equation is given by 



-A 1 







(18) 



which gives A 1.2 = ±-3^ such that the equilibrium point 
(0, 0) is a saddle for A > 0. Since ip is the time derivative 
of cf>, the plot of <f> versus ip will define the phase-plane for 
the equation. The integral curve for the 7^(3, 3) equation 
as obtained from (11) is a hyperbola 
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a constant. 
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Figure 1 gives the phase plane for A = 1 and a particular 
choice of c, say 4. The hyperbolic trajectories reach the 
equilibrium point along two directions only and in all 
other directions the trajectories diverge from it. Thus 
we can say that, in general, the trajectories diverge from 
the critical point making it unstable. 

In Figure 2 we plot the vector (or tangent) field for the 
K(3, 3) equation resulting from the integral curve (19). 
Clearly, at all points the tangent vectors diverge from 
the critical point (0, 0) and reconfirm that it is an un- 
stable equilibrium of the system. For A < 0, the critical 
point becomes a centre. A local bifurcation occurs when 
a parameter change causes the stability of the equilib- 
rium point to change. Thus the system modelled by the 
K(3, 3) equation is expected to exhibit bifurcation as the 
values of A are changed from negative to positive. 

For the critical point at the origin and A > we write 
(0, -0) = ($, Specializing (11) to the K(3, 3) case 
we linearize the system about (0, 0) and write the matrix 
equation 
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FIG. 2: Vector flow for the integral curve in (19) 



Since the eigenvalues of the square matrix in (20) are 

we can write the general solution of the linear equa- 
tion as 
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where ai and a2 are constants of integration determined 
by the initial values $(0) and *(0). The K(2, 2) equa- 
tion in (16a) is similar to that in (16d) except that we 
have an inverted oscillator with shifted origin. It is quite 
straightforward to show that our analysis for the K(3, 3) 
equation also applies to (16a). 

As opposed to K(2, 2) and K(3, 3) equations, there is 
no simple relationship between the K(2, 3) and K(3, 2) 
equations. Curiously enough, reduction of the K(2, 3) 
equation to two first-order equations does not define a 
critical or stable point of the system. Rather, (166) gives 
the motion of a particle with constant acceleration. On 
the other hand, we can easily write such first-order equa- 
tions for the K(3, 2) equation giving the critical points 

(±|y^, 0) for v = 1. An important implication of this 
result is that A should always be less than zero in (16c). 
The choice A = — 1 leads to the critical points (±-^p, 0). 
The eigenvalues of the Jacobian matrix of this problem 
at (^, 0) are purely imaginary while those at (— 0) 
are real and have opposite signs such that the former crit- 
ical point is a centre while the later one is a saddle. The 
integral curve as obtained from the first-order equations 
corresponding to (16c) is given by 



6 3 - -0 + 5i/> 2 = 



(22) 



The solutions tj> — f(ip) of this equation will give the 
phase trajectory. Since (22) is cubic in (p, we can write 
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FIG. 3: Phase diagram for the K(3, 2) equation 



its solutions as 
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FIG. 4: Vector flow for the integral curve in (22) 



which represents a compacton for a% — a 2 — a (say). On 
the other hand, similar linearization of the system about 
the saddle point (—^, 0) leads to a soliton solution. 

4. Conclusion 
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(23c) 

It can be seen that for -0.409682 < ip < 0.409682, all 
(j>iS are real else one of the 0j's is real and other two 
are complex. For < 0.409682 we present, in figure 3, 
the phase trajectories defined by (23). The trajectories 
of a and f3 form a closed curve and clearly show that the 
system has a centre at (^?, 0), while (3 and 7 curves are 
similar to those of figure 1 such that the critical point 
(— 0) is a saddle. Figure 4 portrays the vector field 
corresponding the integral curve (22). Looking closely 
into this figure we see that the tangent vectors converge 
towards the critical point (-^,0) labeled by a and di- 
verge away from (— 0) labeled by b since they repre- 
sent the stable and unstable equilibria of the system. 

Linearization of the system of equations implied by 
(16c) about the centre 0) we get 
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1/3 The nonintegrable K(n, m) equations of Rosenau and 

Hyman have been extensively studied in the literature 
(236) with a view to construct their solutions which appear 
in the form of compacton, peakon, cuspon etc.. In 
this work we did not try to present any new method 
to solve these equations. On the contrary, we made 
use of the phase-plane analysis of dynamical systems 
theory to study the nature of the solutions and their 
interconnection. Thus one may reasonably ask whether 
our study could provide any added realism for the 
solutions of the equations. Our answer to this query 
is fairly straightforward. We worked with K(n, m) 
equations characterized by a parameter A that can 
continuously tune the effects of nonlinearity. We found 
that unlike evolution equations with linear dispersive 
terms the FNE equations do not form Hamiltonian 
systems. However, the later reduce to the Hamiltonian 
form on the zero-energy potential hypersurface. On 
this surface the K(2, 2) and K(3, 3) equations have 
two critical/equilibrium points, namely, a centre and a 
saddle, the presence of which lead to compacton and 
soliton solutions respectively. As the parameter A is 
continuously varied from positive to negative values the 
compacton solution goes over to the soliton solution. 
This represents a signature of a local bifurcation. For 
the K(3, 2) equation the variation of A is constrained 
to negative values only. But as with the other two 
equations the if (3, 2) equation is also characterized by 
two critical points leading once again to compacton and 
soliton solutions. But one cannot be reached from the 
(24) other by varying A. The K(2, 3) equation describes 
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motion of a particle with constant acceleration. Wc 
conclude by noting that these features of the solutions 
did remain undiscovered in any method for solving them, 
however efficient it might be. 
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